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We report a new STAR measurement of the longitudinal double-spin asymmetry All for inclusive 
jet production at mid-rapidity in polarized p + p collisions at a center-of-mass energy of *Js = 
200 GeV. The data, which cover jet transverse momenta 5 < pr < 30 GeV/c, are substantially more 
precise than previous measurements. They provide significant new constraints on the gluon spin 
contribution to the nucleon spin through the comparison to predictions derived from one global fit 
of polarized deep-inelastic scattering measurements. 
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PACS numbers: 21.10, 14.20.Dh, 13.87Ce, 13.88.+e, 14.70. Dj, 13.85.Hd, 12.38.Qk 



The gluon spin contribution to the nucleon spin, AG, 
has been the focus of experimental and theoretical efforts 
since polarized deep- inelastic scattering (DIS) experi- 
ments found that the quark contribution to the nucleon 
spin is small [H . Insight into AG has been obtained from 
next-to-leading-order perturbative quantum chromody- 
namics (NLO pQCD) analyses of the scale dependence of 
the inclusive spin structure function 0, [H, 0, [a, 0, 0] and 
from measurements of hadroproduction of pions and pho- 
tons [1[ . Recent semi- inclusive DIS measurements Q and 
hadroproduction measurements of jets and pions have 
been made 1^, with the latter now being incorpo- 
rated in NLO pQCD fits [12]. Despite recent progress, 
significant uncertainty remains regarding the magnitude 
and sign of AG [HI]. The inclusive measurements pre- 
sented here span more than an order of magnitude in par- 
tonic momentum fraction (x) and are expected to sample 
a sizable piece of the total integral AG = Ag{x)dx. 
Comparisons with predictions derived from one global 
fit [14| to deep-inelastic scattering measurements demon- 
strate the substantial new constraints these results place 
on AG. 

In this paper, we report a new measurement of the lon- 
gitudinal double-spin asymmetry All for mid-rapidity 
inclusive jet production in polarized p + p collisions at 
y^s = 200 GeV center-of-mass energy, 



All = 



(1) 



where a ++ (a^ ) is the differential cross section when 
the beam protons have equal (opposite) helicities. We 
have previously measured the helicity-averaged cross sec- 
tion [lfj for transverse momenta (pr) up to ~50GeV/c 
and it is well described by NLO pQCD evaluations. In- 
clusive jet production in the kinematic regime studied 
here is dominated by gluon-gluon (gg) and quark-gluon 
(qg) scattering. Therefore^-A^L provides direct sensitiv- 
ity to gluon polarization 14] and the cross-section result 
motivates the use of NLO pQCD to interpret our mea- 
surements. 

The data presented here are extracted from an inte- 
grated luminosity of 2pb _1 recorded in the year 2005 
with the STAR detector [l]| at RHIC. The polariza- 
tion was measured independently for each of the two 
counter-rotating proton beams and for each fill using 
Coulomb-Nuclear Interference (CNI) proton-Carbon po- 
larimeters [iH ]. which were calibrated via a polarized 
atomic hydrogen gas-jet target [13] . Averaged over RHIC 
fills, the luminosity weighted polarizations for the two 
beams were 52 ± 3% and 48 ± 3%. The proton helici- 
ties were alternated between successive bunches in one 
beam and between bunch pairs in the other beam. Addi- 
tionally, the helicity configurations of the colliding beam 



bunches were changed between beam fills to minimize 
systematic uncertainties in the All measurement. Seg- 
mented Beam-Beam Counters (BBC) |18J located up and 
downstream of the STAR interaction region (IR) mea- 
sured the helicity dependent relative luminosities, iden- 
tified minimum bias (MB) collisions, and served as local 
polarimeters. 

The STAR subsystems used to measure jets are the 
Time Projection Chamber (TPC) and the Barrel Electro- 
magnetic Calorimeter (BEMC) 0. The TPC provides 
tracking for charged particles in the 0.5 T solenoidal mag- 
netic field for pseudo-rapidities —1.3 < r\ < 1.3 and 27r 
in the azimuthal angle <f>. In 2005 the BEMC, covering a 
fiducial area of </> = 27r and < r\ < 1, provided triggering 
and detection of photons and electrons. 

Events were recorded if they satisfied both the MB 
condition, defined as a coincidence between east and west 
BBCs, and either a jet patch (JP) or high tower (HT) 
trigger. The HT condition required the energy of a single 
calorimeter tower to be at least 2.6 (HT1) or 3.6 (HT2) 
GeV. The JP trigger fired if the sum of a A77 x A(f> = 
lxl patch of towers, the typical size of a jet, exceeded 
4.5 (JP1) or 6.5 (JP2) GeV. Approximately half of the 
2.38 x 10 6 jets extracted from the 12 x 10 6 event set 
originated from the JP2 trigger sample. 

Jets were reconstructed using a mid-point cone al- 
gorithm [l|| with the same parameters as described in 
Ref. [13]. The algorithm clusters TPC charged track mo- 
menta and BEMC tower energies within a cone radius 
of R = A<f> 2 + Ar/ 2 = 0.4. Jets were required to have 
Pt > 5 GeV/c and point between 77 = 0.2 — 0.8 in order to 
minimize the effects of the BEMC acceptance on the jet 
energy scale. BBC timing information was used to select 
events with reconstructed vertex positions within ^60 
cm of the center of the detector, ensuring uniform track- 
ing efficiency and matching the conditions used in de- 
termining the relative luminosity measurements. Beam 
background from upstream sources observed as neutral 
energy deposits in the BEMC were minimized by requir- 
ing the neutral energy fraction of the jet energy (NEF) 
to be less than 0.8. A minimum NEF of 0.1 was also im- 
posed in order to reduce pile-up effects. Finally, only jets 
which contained a trigger tower or pointed to a triggered 
jet patch were considered for analysis. 

Figure Q] compares the NEF spectra for MB, HT, and 
JP triggered jets from data and simulations. Monte Carlo 
events were generated using PYTHIA 6.205 [20j with pa- 
rameters adjusted to CDF 'Tune A' settings 2l| and 
processed through the STAR detector response package 
based on GEANT 3 [H|]. The shapes of the data distri- 
butions are sufficiently reproduced by the simulations for 
the purpose of estimating systematic errors. In contrast 
to the calorimeter triggers, the mean and shape of the 
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MB distribution is relatively stable as a function of jet 
Pt- The HT jets, and to a lesser extent the JP jets, show 
a strong bias towards higher NEF at low px which di- 
minishes for higher jet py. The enhancement of jets near 
NEFwl in the data compared to simulation is consistent 
with contributions from beam background, as discussed 
above. 

We present the inclusive jet All measurement, not as 
a function of the measured transverse momentum (de- 
tector jet pt) : but instead corrected to reflect the jet 
energy scale before interaction with the STAR detector 
(particle jet pt) ■ This correction was carried out by ap- 
plying the same jet reconstruction algorithm to the sim- 
ulated event samples at both the particle and detec- 
tor levels. Particle jets are composed of stable, final- 
state particles which result from the fragmentation and 
hadronization of the scattered partons and remnant pro- 
tons. Detector jets consist of the reconstructed TPC 
tracks and BEMC tower energies in simulated events that 
pass the same trigger conditions that were placed on the 
data. As shown in Fig. [2^a) the jet yield is a rapidly 
falling function of jet pt- This effect combined with the 
jet pt resolution of about 25% results in a shift of the 
inclusive detector jet pt distribution to larger values 
as shown in Fig. Efb). For pt > lOGeV/c the recon- 
structed detector jet pt is on average a factor 1.22 
larger than the PARTICLE jet pt- This shift varies only 
slightly with trigger and underlying partonic process (gg 
vs. qg vs. qq). The dominant uncertainty in the jet pt 
values arises from the ± 4% uncertainty in the jet energy 
scale, but we also account for the subsample dependence 
in the detector to PARTICLE jet conversion. The agree- 
ment between data and simulations is best at high pt (as 
shown in FigfT]) where the correction, proportional to pt, 
is largest. The uncertainty in the particle jet pt slope 
results in a ±2.5% error on the pt scale. 
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FIG. 1: Neutral energy fraction of the jet energy for MB 
(crosses), HT (circles), and JP (squares) data compared with 
STAR simulations for two jet pr bins, (a) 6.2 < pr < 11.4 
GeV/c and (b) 17.3 < pr < 32.2 GeV/c. The statistical 
uncertainties are represented as error bars for the data points 
and bands for the simulations. (Color available on-line) 
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FIG. 2: (a) The raw detected jet yield in data (points) com- 
pared with the STAR Monte Carlo simulations, (b) Corre- 
lation between the reconstructed jet transverse momenta at 
the particle and detector levels. The points indicate the 
means and the vertical error bars show the r.m.s. widths of 
the associated particle jet distributions within the detec- 
tor jet bins. The dashed line represents the condition when 
the particle and detector jet pr values are equal. 



The asymmetry All was evaluated according to, 
E(PiP2)(N++ -RN+-) 



All = 



J2(P 1 P 2 Y(N++ + RN+-)' 



(2) 



in which P\$. are the measured beam polarizations, N ++ 
and denote the inclusive jet yields for equal and 

opposite proton beam helicity configurations, and R is 
the measured relative luminosity. Each sum is over 10 to 
30 minute long runs, a period much shorter than typical 
time variations in critical quantities such as P\,i and R. 
Typical values of R range from 0.85 to 1.2 depending on 
fill and bunch pattern. 

Figure [3] shows the results for inclusive jet All ver- 
sus jet pt corrected for detector response to the particle 
level. The vertical error bars show the statistical uncer- 
tainties. The height of the gray band on each data point 
indicates the total systematic uncertainty on All while 
the width reflects the systematic uncertainty on jet pt- 
An overall 9.4% scale error due to the uncertainty in the 
RHIC CNI polarimeter calibration is not included in the 
systematic error shown in Fig. [3] The present results are 
in good agreement (y /ndf = 7.3/6) with our previous 
measurements of All [101 ] . The combined statistical and 
systematic asymmetry uncertainties are reduced by a fac- 
tor of 4 and the pt coverage is nearly doubled, extending 
up to 30 GeV/c. 

The curves shown in Fig.[3]are NLO pQCD evaluations 
for All based on different polarized parton distribution 
functions [3| • The curve labeled GRSV-std uses the best 
fit to the inclusive DIS data (Ag(x, Qq) = 0.24) at the 
initial scale Q 2 , = 0.4 GeV 2 /c 2 of the parton parametriza- 
tions in Ref. [4|. The other curves correspond to max- 
imally positive (Ag(x,Qo) — g), maximally negative 
(Ag(x, Qq) = —g) and vanishing (Ag(x, Qq) = 0) gluon 
polarization. The calculations were performed with fac- 
torization and renormalization scales [If = [1r = pt- 
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FIG. 3: Longitudinal double-spin asymmetry All for inclu- 
sive jet production at y/s = 200 GeV versus jet pr- The 
points show results for particle jets with statistical error 
bars, while the curves sho w p redictions for NLO parton jets 
from one global analysis [14l ] . The gray boxes indicate the 
systematic uncertainties on the measured All values (ver- 
tical) and in the corrections to the measured jet pr and the 
conversion between particle jet and NLO parton jet pr (hor- 
izontal). (Color available on-line.) 

The A L l values for the GRSV-std and Ag(x,Ql) = g 
cases change by less than 20% for variations in the scale 
frompT/2 to 2px- The calculations are performed for jets 
composed of NLO partons which do not include effects 
due to hadronization and the underlying event. This dif- 
ference, estimated from simulation studies, causes a 1 % 
systematic shift in jet pr between particle and NLO 
parton jets. 

The leading systematic error contribution to the All 
measurement arises from trigger and reconstruction ef- 
fects which cause the asymmetry to differ for particle 
and detector jets. The shift in jet energy scale results 
in the smearing of the particle jet All across the de- 
tected jet pt bin, an effect which is largely accounted 
for by the correction of measured detector jet to par- 
ticle jet px values. The calorimeter triggers, designed 
to select a subset of all minimum bias events, change 
the natural distribution of qq, qg and gg events which 
comprise the inclusive measurement. The consequence 
of this change for All also depends on the true value 
of the gluon helicity distribution. Therefore the system- 
atic error due to both triggering and reconstruction bias 
was estimated from the jet asymmetries calculated within 
the simulation framework for GRSV-std, Ag(x, Qq) = 0, 
and Ag(x, Qq) = — g scenarios. The Ag(x, Qq) — g sce- 
nario, shown in Fig. [31 is not consistent with our data 
and therefore was not included in the estimates. The 
maximum positive and negative differences for the dis- 
tributions were selected at each px bin. Other system- 
atic uncertainties include effects from relative luminosi- 
ties (9 x 10~ 4 ), beam background (7 x 1CV 4 ) and non- 
longitudinal beam polarization components at the STAR 



IR (3 x 1CV 4 ). Parity violating single-spin asymmetries 
in the data were found to be consistent with zero, < 0.2 
standard deviation, as expected, given that parity vio- 
lating physical processes are predicted to be negligible at 
the current level of statistics. 

Figure [4ji illustrates the gluon x range accessed in a 
low and high pr data bin. Note that the sampled Ag(x) 
distributions depend on the magnitude and momentum- 
dependence of the gluon polarization and may have very 
different shapes. The smooth curve represents the cor- 
responding fraction of AG sampled for x m i n = x g i uon in 
the GRSV-std scenario at a scale of Q 2 = 100 GeV 2 /c 2 
which is typical for the present data. The measurements 
presented here provide sensitivity to ^50% of AG. 

To quantify the impact of the new data on AG, the 
measured All values have been compared to predictions 
within the GRSV framework [4|, [lj, [23j in which the po- 
larized parton distributions were re-fit assuming AG is 
constrained to a series of values spanning the full range 
—g < Ag(x, Qq) < g. Figure HJd shows the confidence 
levels (C.L.) found from comparisons with this inclusive 
jet data. The correlations among the systematic uncer- 
tainties for various jet pt have been included in the C.L. 
calculations. 

We find that the presented data exclude fits of AG > 
0.33 at a scale of 0.4 GeV 2 /c 2 with at least 90% C.L. and 
the Ag(x, Q 2 ) = —g scenario is excluded at the 94% level. 
As discussed in Ref. [J] , the GRSV-std fit to the existing 
DIS world data corresponded to a AG{Ql = 0.4) = 0.24 
with a range of —0.45 < AG < 0.7 allowed with a x 2 
variation of 1 . Although these conclusions are dependent 
on the functional form for the gluon polarization defined 
in the GRSV framework, the constraints placed by our 
data on the slice of AG between x=0.02-0.2 are signif- 
icant and will exclude additional PDF's which have a 
AG = J Q Q2 Ag(x)dx contribution larger than GRSV-std 
in this x region. 

In summary, we report new measurements of the lon- 
gitudinal double-spin asymmetry All for inclusive jet 
production at mid-rapidity in polarized p + p collisions 
at = 200 GeV with coverage in jet transverse momen- 
tum up to 30 GeV/c and improved precision compared to 
previous measurements. The data have been compared to 
predictions of All within the GRSV framework, a com- 
monly used set of polarized parton distribution functions, 
constraining AG(Qg) to less than 65% of the proton spin 
with 90% confidence. A global analysis of DIS and RHIC 
data is needed to realize the full impact of these results 
on the shape and magnitude of Ag(x, Q 2 ). 
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FIG. 4: (a) The solid curve represents the fraction of AG for 
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input scale of 0.4 GeV 2 /c 2 4, 14, 23] (Color available on-line). 
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TABLE I: All for inclusive jet production in jet pr bins with 
statistical (stat) and systematic (sys) uncertainties. Both the 
average measured jet pr, (Pt)i an d the final shifted jet pr 
with associated systematics are reported. 



(p¥) 


(p T )+sys.-sys. 


All ± stat +sys-sys 


(GeV/c) 


(GeV/c) 




(xlO- 3 ) 


5.58 


5.60 +0.36 - 


0.41 


-1.51 ± 5.92 +2.11 -2.11 


6.84 


6.14 +0.46 - 


0.48 


-1.16 ± 5.44 +1.77 -1.90 


8.38 


6.83 +0.53 - 


0.54 


3.27 ± 5.59 +1.61 -1.86 


10.26 


8.67 +0.65 - 


0.61 


13.51 ± 6.53 +2.91 -2.15 


12.57 


10.34 +0.83 


-0.76 


-10.01 ± 8.58 +2.60 -1.71 


15.41 


12.89 +0.96 


-0.83 


0.19 ± 12.59 +2.77 -1.78 


18.90 


15.65 +1.09 


-0.88 


-7.30 ± 20.21 +3.09 -2.22 


23.20 


19.30 +1.31 


-1.02 


-39.69 ± 34.84 +2.99 -2.99 


28.39 


23.48 +1.59 


-1.19 


83.17 ± 67.03 +4.10 -4.10 


34.81 


27.94 +1.88 


-1.37 


52.70 ± 137.65 +5.74 -3.73 



